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ABSTRACT 

We present new 0.9-2.45/im spectroscopy (i? ~ 1000), and Y, J, H, Kg, L' photometry, obtained 
at Gemini North, of three low-mass brown dwarf companions on wide orbits around young stars of 
the Upper Scorpius OB association: HIP 78530 B, [PGZ2001] J161031.9-191305 B, and GSC 06214- 
00210 B. We use these data to assess the companions’ spectral type, temperature, surface gravity and 
mass, as well as the ability of the BT-SETTL and Drift-Phoenix atmosphere models to reproduce 
the spectral features of young substellar objects. For completeness, we also analyze the archival 
spectroscopy and photometry of the Upper Scorpius planetary mass companion IRXS J160929.1- 
210524 b. Based on a comparison with model spectra we find that the companions, in the above 
order, have effective temperatures of 2700±100K, 2500±200K, 2300±100K and 1700±100K. These 
temperatures are consistent with our inferred spectral types, respectively M7 (3, M9 7 , M9 7 , and 
L 47 , obtained from spectral indices and comparisons with templates. From bolometric luminosities 
estimated from atmosphere model spectra adjusted to our photometry, and using evolution models at 
5-10 Myr, we estimate masses of 21-25Mjup, 28-70 Mj^p, 14-17Mj^p, and 7-12 Mjup, respectively. 
[PGZ2001] J161031.9-191305 B appears significantly over-luminous for its inferred temperature, which 
explains its higher mass estimate. Synthetic spectra based on the BT-Settl and Drift-Phoenix 
atmosphere models generally offer a good fit to our observed spectra, although our analysis has 
highlighted a few problems. For example, the best fits in the individual near-infrared bands occur at 
different model temperatures. Also, temperature estimates based on a comparison of the broadband 
magnitudes and colors of the companions to synthetic magnitudes from the models are systematically 
lower than the temperature estimates based on a comparison with synthetic spectra. 

Subject headings: brown dwarfs — infrared: planetary systems — stars: imaging — stars: individual 
(HIP 78530, [PGZ2001] J161031.9-191305, GSG 06214-00210, IRXS J160929.1- 
210524) — stars: low-mass 


1 . INTRODUGTION 

Since the first confirmed detection of a brown dwarf 


in 1995 (Oppenheimer et al. 19951, a substantial ef¬ 
fort has been made on characterizing substellar objects. 
Up to now, almost 2000 isolated brown dwarfs have 
been discovered, and we are beginning to better un¬ 
derstand their formation and evolution. Nevertheless, 
the modeling of their cool atmospheres, bearing several 
molecules and dust clouds, is a great challenge for mod¬ 
ern astrophysics. Several low-mass substellar compan¬ 
ions (^30 Mjup) have been discovered recen tly on wide 
orbits (> 80 AU) a round stars, see for example Neuhauser 
& Schmidt| (|2012|) and references therein. Tne age and 
distance of these companions can be inferred from their 
primary star, while their large separation permits their 
direct observation without the hampering glare of their 
primary star; their characterization is thus particularly 
interesting for testing low temperature atmosphere and 
evolution models. At young ages these companions are 
even more interesting as this is where an empirical veri- 
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fication of the theoretical models is most needed (Allard 


et al. 2012). Furthermore, these young companions are 
potentially (more massiv e) analogs to the young giant 
planets recently imaged dMarois et al. 2008 Lafreniere 
et al.||2008[ ILagran^ et al.||2010| IGarson e/al.|^13[ 


Kuzuhara et al.|2013 Rameau et al.|2013p, and thus can 

serve as workbenches in support of the more difficult di¬ 
rect imaging studies of exoplanets. 

In this paper we present and analyze new near-infrared 
photometric and spectroscopic observations of a sam¬ 
ple of three wide substellar companions to young stars 
in the Upper Scorpius OB (USco) formation region: 
the ^ 16 Mjup companion at a separation of ^ 320 AU 
around the K7 star GSG 06214-00210 (hereafter G06214; 
Ireland et al.|2011), the ~ 23 M iuu companion at 740 AU 
trom the H9 star HIP 78530 (Lafreniere et al. |2011 ), 
and the ~ 34 Mjup companion at ~ 885 AU from the 
K7 st ar [PGZ20011 J1610 31.9-191305 (hereafter J1610- 
1913; Kraus et al. 2008). The primary of J1610-1913 
is itself a tight binary (Aab , separation of ^01'145 or 
^20AU, Kraus et al. 2008| ), in which the companion 
(Ab) is roughly at the stellar/substellar boundary. We 
also apply the same analysis, using existing data, to 
the ~ 8 Mjup companion at 330 AU around the K7 star 
IRXS J160929.1 -210524 (hereafte r J160 9-2105) that was 
first ident ified in Lafreniere et al. (2008| and further an¬ 
alyzed in Lafreniere et al. 12010). USco is located at a 
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distance of 145 ± 14 pc (de Zeeuw et al. 1999 fPreibisch 
& Zinnecker 19991 and the average age in the region is 
estimated at 5 Myr, with a very sma ll scatter (±lMyr) 
(De Geus||1992 Preibisch et al.(|2002 1. It is thus reason- 
able to consider a common age tor the systems in our 
sample, meaning that on a comparative basis, the com¬ 
panions should not be affected by the age-mass degener¬ 
acy inherent to substellar objects. The initial age esti- 
m ate of 5 Myr fo r USco was recently revised to 11±2 Myr 
by Pecaut et al. (20121, based on isochrone fitting. This 
new age see ms to be consistent with the results of |Song| 
et al. (|2012 ), given recent the revision of the age of the 
Beta Pictoris moving group proposed by|Binks fc Jeffries] 


Y filter, to complement similar observations made previ¬ 


ously in J, H and Kg and initially reported in Lafreniere 
et ai.|(|2011|). We also obtained observations of HIP 78530 


(2013). Still, the exact age of USco is still the subject of 
debate, and is beyond the scope of this paper, so in our 
work below we present results assuming both estimates. 

The four companions studied in this paper have 
all been studied to various degrees in earlier publica¬ 
tions. The near-infrared photo metry and spectroscopy 
of HIP 78530 B wa s analyzed in Lafreniere et al.|(20Il). 
Bailey et al.| (|2013| studied HIP 78530 B, G06214 B, and 
J1609- 2105 b using 3-5 fim photometry. [Bowler et al.| 
(2011) presented 1.1-1.8/rm spectroscopy of G06214 B 


adding to the JHKL' photometry from Ireland et al. 
(2011). The latter study also independently confirmed 


the common proper motion of J1609-210 5 b. J1610- 
1913 B was observed in the Ks band by |Kraus_e^^T 


( |2008| ), who made its discovery. Recently, [Alier et al. 
(2013) presented a low resolution {R ~ lOOj 0.8-2.5 /rm 


spectrum and H- and iL-band photometric measure¬ 
ments of J1610-1913 B. In this paper, in addition to 
carrying out a homogeneous analysis of the four compan¬ 
ions, with a focus on a comparison of their spectra with 
those of atmosphere models, we also present new data. 
We present the first photometric measurements in Y, J 
and V of JI610-1913 B. For HIP 78530 B, the F-band 
photometry and the 0.90-1.15/rm spectrum have never 
been published before. The spectrum of G06214 B in the 
K band is also presented for the first time, along with the 
part of the J band between 1.00/rm and 1.18 ^im, where 
the VO band and the Na I and the two first K I lines 
are found. We also present the hrst Y -band photometric 
measurement for G06214 B. 

2 . OBSERVATIONS AND DATA REDUGTION 
2.1. Imaging 

The imaging observations were performed at the Gem¬ 
ini North telescope in semester 2011A (program GN- 
2011A-Q-60) using the NIRI camera in co mbination with 
the A LTAIR adaptive optic (AO) system ([Herriot et al. 
2000). The primary stars themselves were used tor wave¬ 
front sensing and the ALTAIR field lens was used to re¬ 
duce the effects of anisoplanatism and achieve better im¬ 
age quality at separations of a few arc seconds. The //32 
camera was used, resulting in a pixel scale of 0.0214" and 
a field of view of 22" x 22" The Cassegrain rotator was 
turned off during the observations, to match the setup 
used for earlier observations of the same stars, and thus 
the field of view orientation changed slowly during the se¬ 
quences. For HIP 78530, we took observations with the 

® As given on the instrument web page at http: //www. gemini. 
edu/sciops/instruments/niri/imaging/pixel-scales-and-fov 


in the narrowband ATcont filter (2.0975 /im) for astromet¬ 
ric follow up as observations at three earlier epochs had 
already been obtained in this filter. For GSC 06214 and 
J1610-1913, we took images with the Y, J, i7, and ATg 
filters. The observation log is presented in table 

For all targets we used a pattern of five dither posi¬ 
tions consisting of the centre and corners of a square of 
10" on one side. For most observations, the primary is 
too bright to get a good signal from the companion with¬ 
out saturating the detector. To obtain deeper images al¬ 
lowing more precise photometry of the companions, we 
thus obtained, at each dither position, a set comprising 
unsaturated images consisting of multiple co-additions of 
short integrations in fast, high read-noise mode, followed 
by one saturated image consisting of one long integra¬ 
tion in slow, low read-noise mode. The saturated images 
can be readily registered to the unsaturated images and 
easily corrected in the saturated part using the properly 
scaled unsaturated images. 

We also observed HIP 78530 and JI610-1913 in L', 
still with the //32 camera but without the use of AO. 
We observed the fai nt photometric standard star FS 140 
( Leggett et al.|2003 ) shortly after the targets to calibrate 
the L' photornetric measurements. For these observa¬ 
tions, we used a pattern of five A — B nod pairs with 
a separation of 8", each pair being displaced from the 
preceding one by 2". At each position we obtained 12 
co-additions with an integration time of 0.75 s each, en¬ 
suring that the primary star was never saturated. This 
sequence was repeated three times for HIP 78530, with 
a pattern rotation of 90° between each sequence. For 
J1610-I913, the sequence was executed only once. 

We reduced the data using custom IDL routines. For 
the images obtained in the high read-noise mode, a 
striped noise pattern was often present and we removed it 
by proper median filtering. For the images taken in F, J, 
H, and A', we constructed a sky frame by taking the me¬ 
dian of the images at all dither positions, after masking 
out the sources in each one. For the images in A', we built 
the sky frame as the mean of the two images obtained at 
the preceding and following dither positions (also after 
masking out the sources). After subtraction of the sky 
frame, we divided the images by a normalized fiat field, 
and we corrected the geometric distortion of the images 
using the prescription given on the instrument webpageQ 
The reduced images at each dither position were then 
registered to place the primary star at their centre, de- 
rotated to a common field orientation, and their median 
was taken. The saturated region of the long-exposure 
images were finally replaced by the properly scaled un¬ 
saturated images. For improved consistency among all 
of our targets, we reprocessed the archival J-, H- a nd 
A'-band data for HIP 78530 B ( Lafreniere et al.||201I ). 


2.2. Spectroscopy 


^ The distortion is given by r' = r + k * r^, where k = 
(1.32 it 0.02) X 10“®, r is the uncorrected distance from the field 
centre in pixels, and r' is the corrected distance from the centre in 
pixels. From http://www.geinini.edu/sciops/instruinents/niri/ 
imaging/pixel-scales-and-fov 
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Table 1 

Observation log 


Target 

Date 

Mode 


Total Integration Time Per Filter (s) 

(^non —saturated! ^saturated) 


HIP 78530 

2011-03-30 

Imaging 

y(25, 50), 

7Ccon.(30, 30) 


J1610-1913 

2011-04-19 

Imaging 

y(i75, -), 

J(75, -), 77(50, -), Ks(ib, 50) 


G06214 

2011-04-19 

Imaging 

y(40, 50), 

J(30, 50), 77(50, -), 77460, -) 


J1610-1913 

2011-06-20 

Spectro 


- (600, -) 


G06214 

2011-06-23 

Spectro 


- (2880, -) 


HIP 78530 

2011-07-03 

Spectro 


- (1350, -) 


HIP 78530 

2011-08-16 

Imaging 



L'(288, -) 

J1610-1913 

2011-08-16 

Imaging 



L'(90,-) 


The spectroscopy observations were made at the Gem¬ 
ini North telescope in the same pro gram as the imag- 
ing, using the GNIRS spectrograph (Elias et al.||2006l 
in cross-dispersed (XD) spectroscopy mode with a t)'.'45- 
wide slit, the 10 lines mm“^ grating and the long blue 
camera with its LXD prism, resulting in a coverage from 
0.885 /im to 2.425 /rm, see observation log in table ^ The 
ALTAIR AO system was also used to improve the spa¬ 
tial resolution and image quality, and thus greatly reduce 
the contamination from the bright nearby primary at the 
position of the companion. Given the wide slit used, the 
spectral resolving power achieved is determined by the 
width of the AO-corrected point-spread function (PSF) 
140-190 mas) and varied between 900 and 1300 de¬ 
pending on target and wavelength. We obtained three 
exposures of 100 to 360 s integration, depending on the 
source, at each of two nod positions along the slit (for 
sky subtraction). W e observed the AO telluric s tan- 
dard star HD 151787 (|Houk & Smith-Moore 19881 im¬ 
mediately after each target to determine and correct for 
the effect of the atmospheric and instrumental transmis¬ 
sions. Observatory standard calibration data (flat field, 
arc lamps) were obtained with each observation. 

We reduced the data using custom IDL routines. First, 
we subtracted the exposures taken at two different nods 
in the slit to remove the majority of the background 
signal, resulting in parallel positive and negative signal 
traces. We then divided the frames by a normalized flat 
field, using a different lamp in the K band than for the 
rest of the spectrum for saturation considerations. A few 
frames also presented a noise pattern of stripes that was 
removed by carefully applying ID iterative median fil¬ 
tering for each quadrant separately. We then rectified 
the traces of each order using cubic interpolation. We 
next corrected a slight spectral shearing by rectifying at¬ 
mospheric lines that were apparent on the frames before 
subtracting the two nod positions. With AO, the shape 
of the PSF is wavelength dependent. We thus fitted 
the trace independently for each spectral pixel along the 
spectrum. We fitted an analytic trace consisting of the 
sum of a Gaussian profile for its core and a Moffat profile 
for the wings. While fitting for the trace, we simulta¬ 
neously fitted (and removed) the potential contribution 
from the primary star and any residual background sig¬ 
nal. The contamination from the primary depends on the 
contrast and separation of the companion and is most im¬ 
portant for G06214 B, with only ^2" separation. The 
flux from each nod position, cleaned from contamination 
and residual background signal, was extracted separately 


using the fitted trace as weight; the flux from the two 
positions were then summed together. The wavelengths 
were calibrated using an Ar arc lamp exposure and the 
different orders of the spectrum were combined by ad¬ 
justing their overlapping sections. Then, we divided the 
target spectra by the total transmission function. The 
latter was determined from the spectrum of the telluric 
standard. Namely, the continuum of the standard star 
spectrum was modelled by a black body function and re¬ 
moved, while its hydrogen absorption lines were fitted by 
a Voigt profile over the appropriate wavelength ranges, 
and then divided out. The median of all individual spec¬ 
tra was taken as the final spectrum and their dispersion 
was used to estimate the uncertainties. Based on the 
achieved PSF FWHM, the effective resolving power in 
the H band for the different spectra were R ^ 1110 for 
HIP 78530 B, i? ~ 930 for G06214 B and i? - 1260 for 
J1610-1913 B. 

3. ANALYSIS AND RESULTS 

3.1. Photometry and astrometry 

For each system, the position of the primary and com¬ 
panion was measured by fitting an elongated 2D gaussian 
function to their PSF. For saturated PSFs, the position 
from the preceding unsaturated frame was used. The 
flux ratio between the companion and primary was calcu¬ 
lated using aperture photometry, with an aperture radius 
set to the radius at which the radial intensity profile of 
the companion falls below the Itr background noise level. 
The contribution of the primary star flux inside the pho¬ 
tometry aperture of the companion was estimated, and 
removed, in the following manner. First, an azimuthally 
symmetric median radial intensity profile of the central 
star was calculated and subtracted from the image. Then 
a similar profile was calculated for the companion in the 
residual image, and this profile was subtracted from the 
original image. This process was then repeated once to 
ensure that the radial profile of the primary was not bi¬ 
ased by the companion. The flux measurement for the 
companion was performed on the original image to which 
we subtracted the modelled flux of the primary star, 
while the measurement for the primary was made on the 
original image minus the modelled flux of the companion. 
These measurements were performed on the combined 
images as well as on the individually reduced frames. 
The uncertainties on the separation, position angle and 
photometry of the companions were determined from the 
scatter of measurements from individual frames. The 
pixel scale, (0'.'0214 pixel“^), was taken from the instru- 
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ment manual and the direction towards North was taken 
from the image headers. By comparing our 2011 mea¬ 
surements with measurements of the same systems made 
at earlier epochs, we noticed that our position angle val¬ 
ues were systematically off by (—0.45 ± 0.04)°, based on 
previous measurements on 3 targets; we thus corrected 
our measurements for this systematic offset and included 
it in the position angle errors. The flux ratios between 
the primaries and companions are given in table along 
with their angular separations and position angles. 

The companion apparent magnitudes were computed 
from the measured contrast ratios in combination with 
the J, H and i^s magnitudes o f the primaries taken from 
2MASS ( Skrutskie et al.|[2006 ). For the J1610-1913 sys¬ 
tem, the photometric nieasurements from 2MASS did 
not resolve the tight binary. We thus corrected the 
photometry of the primary according to our measured 
contrast for the tight binary. The 2MASS magnitudes 
were converted to the MKO system using the tra nsfor¬ 
mation equat ions from the online supplementof |Car- 


penter (2001). The T-band magnitudes of the primaries 


were estimated from template spectra of the appropriate 
spectral type take n from the P ickles Atlas stellar spec¬ 
tral flux librarjpl (|Pickles||l998 ). Using the MKO filter 
profiles and zero points from Tokunaga & Vacca (20051 
and online supplementj^ the atlas spectra were scaled 
to fit our measured fluxes in the JHK^ bands, and then 
integrated over the Y filter to get the synthetic F-band 
flux, and thus the F-band magnitude of the star. For the 
L' band, we measured the magnitudes of the primaries 
and companions directly from our images, using our ob- 
servations of the fai nt photometric standard star FS 140 
([Leggett et al.|200^ for calibration, as mentioned earlier. 
TaHesl3p|4p[5pand present the resulting photometry 
and colors for each system. The colors in these tables 
and in Figure along with the spectra presented below, 
have been corrected for interstellar extinc tion using the 
F JHK -hanA absorption coefficients fro m ICardelli et ah] 
(|1989|) and the A {L)/A{V) ratio f r om |C ox| (|2000[), as- 
summg Ry = 3.1. Carpenter et al.|(|2009|) published ex¬ 
tinction values of Ay = 0.5 tor HIF 78530, Ay = 1.1 for 
J1610-1913 and Ay = 0.0, for J1609-2105. For G06214 


Bailey et al. (20131 found that its extinction is consis- 
tent with Ay = 0.0, which we adopted. We also took 
the extinction into account when measuring the F-band 
photometry using the procedure described above, by fit¬ 
ting extinction-corrected fluxes to the template spectra, 
calculating the synthetic F-band magnitude, and then 
applying the proper correction to our measurement. 

The new ATcont data for HIP 78530 B were used to fur¬ 
ther assess the common proper motion of this compan¬ 
ion with its primary, as s ome doubts about its physical 
association were raised by Bailey et al. ([2013 1, who men¬ 
tioned that it could possibly be an early-M background 
star. At epoch 2011.2422, we measure a separation and 
position angle of 4.527" ± 0.003" and 140.30° ± 0.1°, 
respectively, which are consistent with the values of 


4.529"±0.006" and 140.32°±0.1° measured by Lafreniere 


I ® http: / /www. astro. caltech. edu/$\sini$iiiic/2mass/v3/ 
transfor mations/ | 
I ^ http: //WWW. stsci.edu/hst/observatory/cdbs/pickies7 
atias.html 


et al. (2011) for epoch 2008.3940. Over that time, the 
separation and position angle of a (stationary) back¬ 
ground star would have decreased by 0.034" and 0.82°, 
respectively. Thus our new measurements indicate with 
increased significance IOct) that the companion is co¬ 
moving with the primary. Moreover, the new spectrum 
we have acquired is inconsistent with the compa nion be¬ 
ing an early-M background star (see Sections 3.2 and 


4.11. 


3.2. Spectroscopy 

The newly obtained spectra of HIP 78530 B, G062I4 B, 
and J16I0-I9I3 B are show n in Figure The sp ec- 
trum of JI609-2105 b from Lafreniere et al. (2010) is 
also shown. The average per pixel signal-to-noise (S/N) 
ratios of our three spectra over the whole spectral range 
are ~ 130 for J1610-1913 B, -60 for HIP 78530 B and 
— 30 for G062I4 B. The lower S/N of the latter is due 
to its lower brightness, combined with a relatively more 
important contamination from the primary. The spectra 
display the typical morphologies of young late-M dwarfs, 
with prominent water absorption bands. The spectra 
of the four objects also show a smooth gradation in all 
three spectral bands. From top (HIP 78530 B) to bottom 
(J1609-2I05 b) in the figure, the J-band spectrum shows 
increasingly deeper VO and FeH absorption bands. Fur¬ 
thermore, the slopes of the blue side of both the H and 
K bands become increasingly more pronounced, owing 
to stronger absorption by water vapor. 

We used the method of K Gru z et al. (in prepa¬ 
ration; see [Gruz & Niihez 2007) to assign spectral 


types to the objects presented here. The method con¬ 
sists of a band-per-band visual comparison with field, 
intermediate-gravity and very low-gravity spectroscopic 
templates that were constructed from a median combina¬ 
tion of several spectra that were assigned the same spec¬ 
tral type and gravit y class in the NIR. We verified that 
the clas sifications of Lodieu et al. (2007 ) and Allers fc Liu 
(2013a), based on the H 2 U index from Allers et al. ( j20(J7|) 
and t b ^e H 2 O -I, H 2 O -2 and FeH indices from jslesnick 


|et al.|jl2004 ), generally agreed within one subtype (see 
i'able l^. io summarize, we obtained spectral types of 
M7/3,ld97, M97, and L 47 for HIP 78530 B, J1610- 
1913, G06214 B, and J1609-2105 b, respectively. 

We also ap p lied th e gravity classification scheme of 


http://irtfweb.ifa.hawaii.edu/IRrefdata/iwafdv.html ' 


Allers & Liu (2013a) for moderate-resolution spectra, 
which is based on the strength of the FeH, VO 2 , and 
Hcont spectral indices and the equivalent width of Na I 
and K I lines in the J band at 1.138, 1.169, 1.177 and 
1.253 /im. A score of 0 is given to objects having a value 
within Icr from the mean value of the field dwarf se¬ 
quence, a score of 1 or 2 designates intermediate and 
very low-gravity objects, respectively, where the dividing 
criterion is established to roughly separate objects with 
optical gravity classification of /3 and 7 , and the score is 
replaced by the symbol “n” when the spectrum does not 
cover the spectral range of the index or when the index is 
not appropriate for the spectral type of the object. Ta¬ 
ble presents the equivalent widths and their respective 
gravity score for eac h alkali line. Under the scheme of 
[Allers & Liuj ([2013a[), these four equivalent width scores 
account tor one fourth of the final gravity score, as for 
the two FeH indices. The spectral indices and the final 
gravity class (field, intermediate, or very low gravity) are 























































































Wide substellar companions in USco 


5 


Table 2 

Measured astrometric and photometric parameters 



HIP 78530 B 

J1610-1913 Ab 

J1610-1913 B 

G06214 B 

Angular separation {") 

4.527 ± 0.003 

0.171 ± 0.002 

5.943 ± 0.002 

2.204 ± 0.002 

Position angle (deg)^ 

140.30 ±0.1 

90.6 ±0.4 

113.77 ± 0.08 

175.97 ± 0.05 

Ay (mag) 

9.5 ±0.3 

2.78 ± 0.03 

4.64 ± 0.03 

7.00 ±0.26 

AJ (mag) 

8.28 ±0.05'° 

2.54 ± 0.06 

4.02 ±0.02 

6.18 ±0.03 

AH (mag) 

7.61 ± 0.03'° 

2.45 ±0.02 

4.11 ± 0.02 

6.19 ±0.02 

AKs (mag) 

7.28 ±0.03*° 

2.51 ± 0.06 

3.85 ±0.02 

5.74 ± 0.01 

Aiy2_^0? (mag) 

7.27 ± 0.07 




AL' (mag) 

6.9 ± 0.2 

2.50 ± 0.05 

3.33 ±0.04 

4.75 ±0.05° 


a 

b 

c 


Corrected for a —0.45 it 0.04° offF jet 
Rem easured from observa tions of 
From [Ireland et al.| ||2011|| . 


Lafreniere et al. 


( 20111 . 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength (jum) 


Figure 1. In black from top to bottom, our GNIRS-XD spectra of HIP 78530 B, J1610-1913 B, and G06214 B, and the archival spectrum 
of J1609-2105 b from|Lafreniere et al.|(|2010||. The spectra of HIP 78530 B and J1610-1913 B have been corrected for extinction (see text). 
Regions of strong te lluric absorption nave been greyed out. In green from top to bottom, comparison spectra of USco brown dwarfs from 
[Lodieu et ^ p007| : USco J155419-213543 (MS), USco J160830-233511 (M9), USco J160714-232101 (LO), and USco J163919-253409 


presented in TableJ^ For HIP 78530 B, we obtain a score 
of ln21 for the FeIvVOz, Alkali lines and iJ-cont indices, 
respectively, which classifies it as an intermediate-gravity 
object. J1610-1913 B and G06214 B are classified as very 
low-gravity objects, with both a scores of 2n22, respec¬ 
tively. Finally, J1609-2105 b is also classified as very 
low gravity, but with some reserve considering its score 
of 2nnn, the spectral range and low S/N preventing us 


from using all indices but FeHj. The index-based gravity 
classes of all objects agree with our visual classification. 



throughout this work for simplicity. 
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Table 3 

Properties of HIP 78530 


Value 


Parameter 

Primary 

Companion 

Y (mag) 

6.766 ± 0.020^ 

16.27 ± 0.05 

J (mag) 

6.925 ± 0.021*^ 

15.21 ± 0.05 

H (mag) 

6.931 ± 0.029*^ 

14.55 ±0.04 

Ka (mag) 

6.900 ± 0.020*^ 

14.18 ± 0.04 

L' (mag) 

6.91 ± 0.02 

13.81 ± 0.20 

J - Ks (mag)*^ 

-0.57 ± 0.03 

0.95 ± 0.06 

H - Ka (mag)‘= 

0.00 ±0.04 

0.34 ± 0.06 

Ka - L' (mag)^ 

-0.04 ±0.03 

0.3 ±0.2 

Spectral type 


M7±0.5/3 

Teff (K) 

~ 10 500® 

2700±100 

Distance (pc) 

156.7 ± 13.0f 


Projected separation (AU) 

740 ± 60 

log(L/L©) 


-2.53 ±0.09 

Mass (Mq) [5 Myr] 

~2.5® 

0.022 ± 0.001 

Mass {Mq) [10 Myr] 

~2.5® 

0.023 ±0.002 


the measured flux in other bands, see 
^ From 2MASS PSC llSkruts kie et ah. 
system with the equations m|Uarpenter|(|z 001 ||. 
^ Dere ddened colors, see text tor det ail. 

^ FromiHouk Smith-Moo^ 


200 6 I, con verted to the MKO 


® From patremere et al.| ||2011 1 and references therein. 
^ From [^n Leeuwen| ( |2007| ^ 


119881 


Table 5 

Properties of GSC 06214-00210 


Parameter 

Value 

Primary 

Companion 

V (mag) 

10.20 ± 0 . 020 ‘‘ 

17.20 ± 0.05 

J (mag) 

9.946 ± 0.027’^ 

16.13 ±0.04 

H (mag) 

9.329 ± 0.024*" 

15.52 ±0.03 

Ka (mag) 

9.129 ±0.021*> 

14.87 ± 0.02 

L' (mag) 

9.10 ±0.05 

13.75 ±0.07 

J - Ka (mag)® 

0.72 ±0.03 

1.16 ±0.04 

H - Ks (mag)® 

0.63 ±0.04 

0.61 ± 0.04 

Ka - L' (mag)® 

-0.01 ± 0.05 

1.08 ±0.07 

Spectral type 

K7±0.5'* 

M9±0.5 7 

Teff (K) 

4200 ± 150'* 

2300±100 

Distance (pc) 

145 ± 14® 


Projected separation (AU) 

320 ± 30 

logfL/L©) 

-0.42 ±0.08'* 

-3.01 ± 0.09 

Mass {Mq) [5 Myr] 

0.9 ±0.1'* 

0.015 ± 0.001 

Mass (Mq) [10 Myr] 

0.9 ±0.1'* 

0.016 ± 0.001 


^ Extrapolated from a template spectrum (|Pickles|1998|l scaled to 
the measured flux in o ther bands, see text to r detail. 

^ From 2MASS PSC jSkrutskie et al.||2006|l, converted in MKO 
with equations in|uarpenter| (|2UUl||. 

^ Dereddened colors, see text tor del 
^ From 


Bowler et ah 
distance ot L 


® Mean _ 

tainties discussed in |Ireland et ai.| < | 2 Ull| |. 


"d^etaih 

l| 2011 l. _ _ 

SCO tromlde Zeeuw et all] )|1999ll, with uncer- 


Table 4 

Properties of [PGZ2001] J161031.9-191305 


Table 6 

Properties of IRXS J160929.1-210525 


Parameter 


Value 


Primary 

Secondary 

iertiary 

Y (mag) 

10.274 ± 0.020® 

12.65 ± 0.05 

14.73 ±0.05 

J (mag) 

10.062 ± 0.026*" 

12.61 ± 0.05 

14.09 ±0.05 

H (mag) 

9.337 ± 0.022*" 

11.80 ±0.04 

13.43 ±0.04 

Ka (mag) 

9.068 ± 0.021*" 

11.58 ±0.04 

12.92 ±0.04 

L' (mag) 

8.72 ±0.07 

11.22 ± 0.07 

12.05 ±0.06 

J - Ks (mag) ® 

0.81 ± 0.03 

0.84 ± 0.06 

0.98 ±0.06 

H - Ka (mag) ® 

0.20 ±0.03 

0.15 ± 0.06 

0.47 ± 0.06 

Ka - L' (mag) ® 

0.27 ± 0.07 

0.29 ±0.08 

0.80 ±0.07 

Spectral type 

K7 '* 

M4 

M9±0.5 7 

Teff (K) 

~ 4000 

3200±300 

2500±200 

Distance (pc) 

145 ± 14 ® 



Proj. sep. (AU) 


26 ±3 

885 ± 85 

log(L/L 0 ) 


-1.48 ±0.11 

-2.13 ±0.12 

Mass (Mq) [5 Myr] 

~0.77 

0.12 ± 0.02 

0.032 ± 0.004 

Mass {Mq) [10 Myr] 

~0.77 ** 

0.16 ± 0.02 

0.058 ±0.011 



Value 

Parameter 

Primary 

Companion 

J (mag) 

9.764 ± 0.027® 

17.85 ±0.12*" 

H (mag) 

9.109 ±0.023® 

16.86 ±0.07*" 

Ka (mag) 

8.891 ± 0.021® 

16.15 ±0.05*" 

L' (mag)® 

8.73 ±0.05 

14.8 ±0.3 

J - Ka (mag) 

0.87 ± 0.03 

1.7± 0.1 

H - Ks (mag) 

0.22 ±0.04 

0.71 ± 0.09 

Ka — L' (mag) 

0.16 ±0.05 

1.4 ± 0.3 

Spectral type 

K7 V±l® 

L4dil 7 

Teff (K) 

4060t“Oc 

1700 ± 100 

Distance (pc) 

145 ± 14'* 


Projected separation (AU) 


330 

log(L/L 0 ) 

-0.37 ± 0.15® 

-3.5 ±0.2 

Mass (Mq) [5 Myr] 

0.85tO;?°® 

0.008 ± 0.001 

Mass {Mq) [10 Myr] 

0.85t°;fo® 

0.011 ± 0.001 


^ Extrapolated from a template spectrum l|Pickles|1998|l scaled to 
the measured flux in other bands, see text for detail. 

^ Resolved MKO photometry based on our measured contrast and 
unr^olved 2MASS PSC photometry jSkrutsk ie et al'1|2006ll, ta - 
ble[ 2 l using the system conversion equations m|Carpenter|(| 2 U 01 ||. 
^ Dere ddened colors, see text for detail. 


'* FromiPreibisch et al. I 20011. 

Mean distance ot 1 
tainties discussed in 

JSco tromlde Zeeuw et ahNl999L with uncer 

Ireland et al.|||L^U11||. 

^ From Kraus et ah 

[T 0081 . 


The Paschen-/? line at 1.282 is detected in emission 
in the spectrum of G06214 B (see Figure]^. Th is emis¬ 
sion line was previously observed and discussed in |Bowl^ 


et ah (2011). Bowler et ah conclude that this emission 
IS a sign of accretion or outflow, revealing the presence 
of a circumplanetary disk. The presence of a disk is con- 


^ From 2MASS PSC llSkrutskie et ah|2006ll, converte d to the 
MKO system with the equations m|<Jarpenter| || 2001 ||. 


^ Base d on the contrasts given in 
From [Lafreniere et ah| | |2008[ |. 
Mean distance of USco from 


Lafreniere et ah 


12008 k 


de Zeeuw et ah 


uncertainties discussed in [Ireland et al.| | | 2 Uil| |. 


119991, with 


sistent with the K — L' excess (1.18 ± 0.10) observed by 


Ireland et al. (20111. Bowler et al. reported an eqnivalent 
width of —11.4±0.3 A, which is significantly higher than 
the value that we measure here (—4.4±0.5 A). This might 
be a sign that the accretion or outflow is variable. We 
verified that this feature is present in our raw spectrum 
(Figure]^, rather that an artifact that could have been 
introduced by an improper correction of the Paschen-/? 
line in our AO telluric standard star. The Brackett-y line 
at 2.166/rm is also present in emission in our spectrum 
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Table 7 
Spectral type 


Object 

SpT 

Visual 


SpT (index) 


H20=‘ 

H 2 O-U 

H20-2‘> 

FeH*^ 

HIP 78530 B 

M7d 

M7.4 ± 0.5 

M8.9 ± 1.2 

M7.5 ± 0.5 

M6.5 ± 1.5 

J1610-1913 B 

M9 7 

M8.4 ± 0.4 

M9.5 ± 1.1 

M8.6 ± 0.5 

M8.3 ± 1.5 

G06214 B 

M9 7 

M9.8 ± 0.6 

L0.2 ± 1.2 

M9.5 ± 0.6 

M9.8 ± 1.5 

J1609-2105 b 

L4 7 

L2.9 ± 1.0 

L2.6 ± 1.2 




a 

b 


From 

From 


^ler^^^J_|200^ 


Slesnick et al. {2004 \ 


Table 8 

Gravity scores from equivalent width 


Object SpT EW (A) [Gravity Score^] 

Na I 1.138 fim K I 1.169 /^m K I 1.177 //m K I 1.253 fim 


HIP 78530 B 

M7 

7.5 ± 0.9 

[1] 

1.2 ± 0.8 

[2] 

2.5 ± 0.7 

[2] 

2.2 ± 0.6 

[1 

J1610-1913 B 

M9 

5.2 ± 0.6 

[2] 

0.6 ± 0.3 

[2] 

1.0 ± 0.3 

[2] 

0.9 ± 0.1 

[2 

G06214 B 

M9 

8.7 ± 1.2 

[1] 

2.5 ± 1.5 

[2] 

5.1 ± 1.3 

[1] 

1.6 ± 1.2 

[2 

J1609-2105 b 

L4 

... 1 

nl 

... I 

nl 

... 1 

nl 

... 1 

■ n 


^ See |Allers fc Liu| ( |2013a| | 


Table 9 
Gravity class 


Object SpT Index Values [Gravity Score^] Gravity 

FeH^ FeHj VO^ Tf-cont Score^ Class Visual 


HIP 78530 B 

M7 

1.062 ± 0.001 

1 


n 

1.072 ± 0.002 

n 

0.987 ± 0.002 

1 

ln21 

INT-G 

J1610-1913 B 

M9 

1.053 ± 0.001 

2 

1.053 ± 0.008 

2 

1.115 ± 0.002 

n 

1.019 ± 0.001 

2 

2n22 

VL-G 

G06214 B 

M9 

1.136 ± 0.005 

1 

1.09 ± 0.03 

2 

1.232 ± 0.005 

n 

1.017 ± 0.001 

2 

2n22 

VL-G 

J1609-2105 b 

L4 


n 

1.04 ± 0.13 

2 


n 


n 

2nnn 

VL-G 


^ See jAllers Sz Liu| ( |2013a| | 

^ Respectively the scores for FeH (highest of FeH^ and FeHj), VOz, alkali line (rounded mean of Na I and K I line scores, see 
Table [^, and 77-cont. 



Figure 2. Zoom on the Paschen-/3 line (1.282 ^m) in our spec¬ 
trum of G06214 B. We measure an equivalent width of (—4.4 ± 
0.5) A. The red error-bar represents the average noise level in the 
plotted region and the green spectrum is the data before telluric 
correction. 


of G06214 B (EW= —0.24 ± 0.05 A), providing further 
evidence for the presence of accretion or outflow. 


4. COMPARISON WITH ATMOSPHERE MODELS 

In the following sections, we compare the spectra of 
the wide companion s in USco to the s ynthetic spectra 
from the B T-Settl (Allard et al.|2011 ) and the Drift- 
Phoenix (Helling etal. 20081 Witte et al. 12009 |20I1 1 
models. Synthetic spectra with Jeff ranging from 1500 K 
to 3500 K, logg ranging from 2.5 to 6.0, and solar metal- 
licity were considered for the analysis and were binned 
to the same spectral resolution as our observed spectra. 
A first fit was performed by minimizing the goodn ess-of- 
fit statistic (G/j) described in|Cushing et al.|(2008D. The 
minimization was performed 10 000 times, each tirrie with 
a Gaussian distribution of random noise, corresponding 
to the uncertainties of our measured fluxes, added to 
our data in each resolution element. The fraction of the 
Monte Carlo simulations (/mc) in which the synthetic 
spectrum was identified as the best fitting model is then 
considered to evaluate the precision of the determination 
of Teff and logg. The fit is evaluated for the whole spec¬ 
trum at once, but also for each band separately. The 
best f its found through this method are summarized in 
Table In addition, the same sets of synthetic spectra 
were visually compared with our observations. The use 
of solar metallicity models for USco is reasonable in light 
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of the results of Mohanty et al. (2004), in particular see 
their section 4.331 We have nevertheless included BT- 
Settl models at higher metallicity {[M/H\ = +0.5) and 
verified that our main conclusions about the companions’ 
properties remained valid. 

We also compare our photometric measurements to 
synthetic magnitudes calculated from the two sets of 
synthetic spectra. To compute the synthetic magni¬ 
tudes, we used t he filter profile and the magnitude 
zero points from [Tokunaga fc Vac^ (20051 and online 
supplementj^ We compared the synthetic and observed 
magnitudes and determined the best-fit model by min¬ 
imizing the over the YJHKsL' bands. In Figure 
we also compare the observed and model colors. As visi¬ 
ble on the figure, the measured colors of the companions 
agree reasonably well with the colors expected from the 
models. In particular, the companions colors seem to 
roughly reproduce the shapes of the model curves. The 
only noteworthy discrepancy is a systematic offset of up 
to ~ 0.1 mag in J—H. The relative positions of the colors 
of the companions in the different color-color diagrams, 
when compared to the model curves, readily indicate the 
relative temperatures of the companions. From the two 
rightmost columns of the figure, we get, respectively from 
the hottest to the coldest, HIP 78530 B, J1610-1913 B, 
G06214 B and J1609-2105 b. This ordering is consis¬ 
tent with the spectral types presented above. In almost 
all panels of Figure we can also see that the colors of 
the companions are closer to the models of low surface 
gravity, as expected for young objects. The only panel 
where this is not the case is J — H vs H — Kg, although it 
seems that this problem would disappear if the ~ 0.1 mag 
systematic offset in J — H mentioned earlier could be re¬ 
solved. 

The temperature estimates based on all of these anal¬ 
yses, along with comments on th e agreement with the 
models, are discussed in Sections |4.l| |4.2[ |4.3[ and |4.4| 
for each object separately. 

4.1. HIP 78530 B 

Figure |5 compares the observed spectrum of 
HIP 78530 B to synthetic spectra selected from the grid 
of BT-Settl and Drift-Phoenix models. For the BT- 
Settl models, the shape of both the J and K bands are 
better fitted by the 2600-2700 K models, at a low grav¬ 
ity (logg = 3.0) to match the CO lines depths, but the 
H band is then too triangular. An effective tempera¬ 
ture of 2800 K is needed to get the right slopes in H 
band. The Drift-Phoenix models at 2600 K and low 
gravity are able to better reproduce the features in the 
J band, particularly the VO band at 1.06/xm. At higher 
temperatures (2800K), this feature is not deep enough 
in the models. The K band is well htted in the 2600- 
2800 K temperature range, with a marginally better fit at 
2600 K. For the H band, a temperature even higher than 
2800 K is needed to correctly ht the shape of the pseudo¬ 
continuum with the Drift-Phoenix models. Using the 
goodness-of-fit statistics (Table 10), the best fits for the 
individual bands are for 2600-29D0K, in good agreement 
with the above estimates, while it is 2600 K for the ht to 
the entire JHK spectrum. As for the broadband pho- 


^ http://irtfweb.ifa.hawaii.edu/'nsfcam/filters.html 
http://irtfweb.ifa.hawaii.edu/IRrefdata/iwafdv.html 


tometry only, the best ht is achieved with of 2300- 
2700 K and log g = 3.5. Considering all of these elements, 
we assign a temperature of 2700 + 100 K to HIP 78530 B. 

The best hts discussed above occur for models at low 
surface gravity, in agreement with the young age of the 
region and with the values of gravity-sensitive spectral 
indices found earlier. In particular, the spectral indices 
for the FeH molecular bands at 0.998/xm and 1.200/xm 
are signihcantly weaker than those of held dwarfs, a sign 
of low surface gravity. As mentioned previously, the VO 
band at 1.06 /xm is gravity-sensitive. Systematically, the 
Drift-Phoenix models provide a much better ht of this 
feature than the BT-Settl model, although they do not 
signihcantly discriminate the surface gravity parameter. 
The depth of the CO molecular bands in the red part of 
the K band does require a low gravity to be well htted. 
The gravity-sensitive Na 1 doublet at 2.206 and 2.209 A 
is clearly visible in the data, but a logg of 4.5 is not 
high enough for the models to reproduce its depth. The 
models are thus under predicting the depth on the Na 
feature. _ 

Based on the photometry of HIP 78530 B, [Bailey et al.| 
(|2013|) estimated a temperature of 3300K-3400 K and a 
spectral type of ^ M3, raising the possibility that the 
companion was instead a background star. While this 
higher temperature estimate could be consistent with the 
K — L' color of the companion, according to models, it 
would not be appropriate for the other colors (see Fig¬ 
ure 1^ . Also, this estimate is inconsist ent with our ob¬ 
served spectrum and the spectrum from Lafreniere et al.| 
( 2011| ). In addition, our observations show signs of low 
gravity, common proper motion (see 3.11, and a spectral 
type of M7. For all of these reasons, we rule out the 
possibility that this is a background star. 


4.2. [PGZ2001] J161031.9-191305 B and Ab 

A comparison of the spectrum of J1610-1913 B with 
various model spectra is shown in Figure]^ The K band 
is best fitted by the BT-Settl model at 2300 K, at low 
surface gravity (log g = 3.0), and by the Drift-Phoenix 
model at 2500 K, for either surface gravities shown. Hot¬ 
ter models (2700 K) fail to match the blue side of this 
band. In the H band, the best fit occurs for temper¬ 
atures of 2500-2700 K, with a notably better ht at low 
surface gravity for 2700 K. In the J band, the best ht of 
the pseudo-continuum as well as the depth of the water 
band at 1.33/xm seems to take place at 2500 K, for ei¬ 
ther surface gravities shown, although the VO band at 
1.06 /xm is not quite deep enough at this temperature in 
the models. The depth of the VO feature is best matched 
by the 2300 K Drift-Phoenix models, for either surface 
gravities, but then the ht is not as good in the other parts 
of the J band. Again, the BT-Settl models fail to re¬ 
produce the VO band at any temperature or gravity. The 
best hts for the individual bands using the goodness-of-ht 
method (Table [l0|) are found for temperatures of 2600- 
2700 K, and the nt to the entire JHK spectrum indicates 
a temperature of 2300-2400 K. All of these best hts oc¬ 
cur for a logg equal to or less than 4.5. For the ht of the 
broadband magnitudes, after proper correction for ex¬ 
tinction as mentioned previously, we get a temperature 
of 2300 K with a log g of 3.0 for both models; this is close 
to the simultaneous JHK model ht result. Considering 
all of these values, we assign an effective temperature 
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Table 10 

Best fit models based on the goodness-of-fit statistics 


Band 

ModeP 

HIP 78530 B 

J1610-1913 B 

G06214 B 


J1609-2105 b 



T’eff 

logg 

/mc 

Ted 

logg 

/mc 

TeS 

logg 

/mc 

'Gff 

logg 

/mc 

j 

BT-S 

2800 

3.5 

0.93 

2600 

3.5 

1.00 

2500*= 

3.5*= 

0.94 

1600 

3.5 

1.00 

j 

D-P 

2900 

5.0 

1.00 

2700 

4.5 

1.00 

2500 

5.0 

0.83 

1600 

3.0 

1.00 

H 

BT-S 

2900 

4.5 

0.61 

2700 

3.5 

1.00 

2600 

3.5 

0.50 

1600 

4.0 

1.00 

H 

D-P 

2900 

4.5 

1.00 

2600 

3.5 

1.00 

2600'= 

3.0== 

0.69 

1800 

3.0 

1.00 

K 

BT-S 

2600 

2.5 

0.88 

2700 

3.5 

1.00 

2600 

2.5 

0.94 

1700 

4.0 

1.00 

K 

D-P 

2800 

4.5 

0.67 

2600 

4.0 

0.78 

2300 

3.0 

1.00 

1800 

3.5 

1.00 

JHK 

BT-S 

2600 

3.0 

1.00 

2400 

3.0 

1.00 

2100*= 

3.0 

1.00 

1600 

3.5 

1.00 

JHK 

D-P 

2600 

3.5 

1.00 

2300 

3.0 

1.00 

2100 

3.0 

1.00 

1600 

3.0 

1.00 

^ BT-S: BT-Settl, D-P:Drift-Phoenix 
^ The VO band was omitted from the fit. 

^ There is also a local minimum at 1700 K but this model is 

clearly not appropriate for other bands, we 


have thus restricted the range of temperatures for the fit to > 1800 K. 



4.5 (green), and 5.5 (orange), and for temperatures ranging from 3000K to 1600K by lOOK increments. The solid lines on the bottom 
row shows synthetic colors from the Drift-Phoenix models for logg of 3.0 (blue), 4.0 (red), 5.0 (green), and 6.0 (orange), for the same 
temperatures. 


of 2500 ± 200 K to J1610-1913 B. The spectral and the 
broadband photometry fits both favour a very low sur¬ 
face gravity for this object, in good agreement with the 
values of the FeH index and alkali lines equivalent widths 
calculated earlier. 

We have not observed the closer-in companion (Ab) in 
the J1610-1913 system using spectroscopy, as getting a 
contamination-free spectrum of this object with a source 
^3 mag brighter at a separation of only ^0.2" is too 
challenging for the instrumental setup we used. Nev¬ 


ertheless, we have photometric measurements from our 
imaging and we can compare those with the models to 
assess its effective temperature. With the photometric 
points obtained, the best fit to the synthetic magnitudes 
would indicate an effective temperature of 2900 — 3300 K 
with a logg of 3.5. 

4.3. GSC 06214-00210 B 

Figurej^shows our GNIRS spectrum of G06214 B com¬ 
pared with a selection of models from the BT-Settl 
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Wavelength (urn) Wavelength (wm) 


Wavelength (urn) 


Figure 4. GNIRS spectrum of HIP 78530 B (black) corrected for an extinction of Ay = 0.5 and compared with synthetic spectra of 
BT-Settl at logg of 3.0 (red) and 4.5 (blue) and Drift-Phoenix at logg of 3.0 (yellow) and 4.5 (green) for Tefr=2600K (top row), 
Teff=2700K (middle row), and Teff=2800K (bottom row). 
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Figure 5. GNIRS spectrum of J1610-1913 B (black) corrected for an extinction of Ay = 1.1 and compared with BT-Settl synthetic 
spectra for logg of 3.0 (red) and 4.5 (blue) and Drift-Phoenix synthetic spectra for logg of 3.0 (yellow) and 4.0 (green), for Teff=2300K 
(top row), Teff=2500K (middle row), and Teff=2700K (bottom row). 
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Wavelength (urn) Wavelength (wm) Wavelength (urn) 


Figure 6. GNIRS spectrum of G06214 B (black) compared with BT-Settl synthetic spectra for logg of 3.0 (red) and 4.5 (blue) and 
Drift-Phoenix synthetic spectra for logg of 3.0 (yellow) and 4.5 (green), for TefF=2200K (top row), Teff=2500K (middle row), Teff=2700K 
(bottom row). 
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and Drift-Phoenix models. The K band is best repro¬ 
duced by models at 2200 K and low logg for both BT- 
Settl and Drift-Phoenix. For the H bands, mod¬ 
els of 2500-2700 K provide reasonable fits, with little ef¬ 
fects from surface gravity. In the J band, the VO and 
FeH (1.2 /rm) bands are better matched by the models at 
2200 K, with Drift-Phoenix providing a much better 
fit than BT-Settl. The most important difference is the 
deeper water absorption band at 1.33/rm for the models, 
as compared with the observations. The goodness-of-fit 
evaluation (Table [To| indicates temperatures of 2300- 
2600 K when appliea to individual bands, and 2100 K 
when applied globally. The colors of this companion 
and the corresponding photometric magnitudes are in 
excellent agreement with the models for a temperature 
of 2200K and very low logg. We assign a temperature 
of 2300 ± 200 K to G06214 B. 

Visually, the models with logg of 3.0 or less are in 
better agreement with the observed spectrum, especially 
in the K band. The same result is obtained for the fit 
of the spectrum and the broad band fluxes. Also for this 
object, the depth of the gravity-sensitive K I doublet at 
1.244 and 1.252/rm in the models is insufficient to agree 
with the observations, even if properly degraded to the 
resolution of the observations. 

4.4. IRXS J160929.1-210524 h 

Fig ure [t] presents the NIFS J band dLafreniere et al . 
2010[ ) anoTMIRI H and K bands (|Latreniere et ai.||20U8|) 
spectrum of J1609-2105 b in comparison to synthetic 
spectra with Tgg ranging from 1600 K to 2000 K, from 
both the BT-Settl and the Drift-Phoenix models. 
The lower-gravity (logg = 4.0) Drift-Phoenix model 
at 1800 K gives the overall best fit, although the red 
side of the J band is a bit too low and the slope on 
the red side of the H band is a bit too steep. The 
1800 K BT-Settl model does not provide as good a fit 
as the 1800 K Drift-Phoenix model, particularly in the 
H band where it is too peaked compared with the ob¬ 
served spectrum. At lower temperatures (Teff=1600K), 
both models clearly fail to reproduce the observations in 
all bands. At higher temperatures (Teff=1800K), the fits 
are not too bad for both models, although the water ab¬ 
sorption band in J is too strong in the models and the 
red side of the H band is too steep. The goodness-of- 
fit evaluations (Table [l0|), both band-by-band and over 
the whole spectrum, generally agree on a temperature 
of 1600-1800K and a logg of 3.0 — 4.0. Only the H 
and K bands fit a higher temperature of 1800 K with 
the Drift-Phoenix models. As for the broad band 
magnitudes, they yield a best-fit temperature of 1700 K 
for BT-Settl and 1800 K for Drift-Phoenix, in both 
cases with a logg of 3.0. We thus assign a tempera¬ 
ture of 1700 ± 100 K to J1609-2105 b. The spectral in¬ 
dices calculated earlier classified this object as having a 
very low gravity; this is in good agreement with the best 
fits with the models obtained here. However here again, 
the gravity-sensitive potassium lines in J are not deep 
enough for the Drift-Phoenix model at lower gravity. 

4.5. Mass estimates 

We have estimated the masses of the companions us¬ 
ing two different approaches, comparing either their es¬ 


timated bolometric luminosities or their estimated effec¬ 
tive temperatures to the predictions of evolution models. 
We have used two s e ts of evolution models, the models 
from Burrow s et al. (19971 and the models from Allard 
(|2U13|), the latter being based on the CIFISi‘2Ull 


et al. 


BT'-Settl atmosphere model s and the AMES-Cond 
isochrones (Baraffe et al. 2003). For the effective tem¬ 
peratures, we simply used the estimates presented in the 
previous section. For each object, a synthetic bolometric 
luminosity was computed for every synthetic spectrum 
within the range of plausible temperature and gravity 
determined previously. This was done by first scaling 
the model spectrum to the observed flux level, by min¬ 
imizing the between the synthetic magnitudes of the 
model spectrum and the ones we observed, and then by 
integrating the entire model spectrum and converting the 
total fl ux to luminosity using the mean USco distance of 
145 pc (|de Zeeuw et al.||1999|) and 14 pc uncertainty (as 
discussed in Ireland et ai.|2(Jll ). The error on luminosity 
is chosen to be large enough to include all the tempera¬ 
tures and gravity ranges described above and encompass 
results from both atmosphere models. The resulting lu¬ 
minosities are included in Tables iii and[^ see also 
Figure for a comparison of these luminosities with the 
prediction of evolution models. 

Estimating the masses based on the calculated lumi¬ 
nosities and an age of 5 Myr, we obtain 0.008 ±0.001 Mq 
for J1609-2105 b, 0.015 ± 0.001 for G06214 B, 
0.032± 0.004 M© for J1610-1913 B, 0.022± 0.001 M© for 
HIP 78530 B, and 0.12 ± 0.02 M© for J1610 Ab. Eor an 
age of 10 Myr, the corresponding masses are respectively, 
0.011 ± 0.001 M©, 0.016 ± 0.001 M©, 0.058 ± 0.011 M©, 
0.023±0.002 Mq, and 0.16± 0.02 Mfr^. The impac t of the 
revised older age proposed by Pecaut et al. (2012) is rela¬ 
tively small for the three lowest mass objects, G06214 B 
and HIP78530 B being in a relatively stable deuterium 
burning phase. The impact for J1610-1913 Ab and B 
would be more important. 

With the objects ordered according to increasing effec¬ 
tive temperature, as above, it is obvious that something 
is off for J1610-1913 B. Namely, its luminosity is much 
higher than expected. Indeed, J1610-I913 B has an esti¬ 
mated temperature of 2500 K, cooler than HIP 78530 B 
at 2700K, but its luminosity (log(L/L©) = —2.13) 
comes out significantly brighter than that of HIP 78530 B 
(log(L/L©) = —2.53). A possible reason for this dis¬ 
crepancy is that the true (unknown) distance of J1610- 
1913 B differs largely from the mean distance of USco 
members. If its true distance were toward the closer side 
of the association, at ~ 115 pc, rather than the assumed 
distance of 145 pc, then its luminosity would be closer 
to log(L/L©) = —2.45. Another possibility is that the 
wide companion is itself an unresolved equal-mass bi¬ 
nary. That would bring its intrinsic luminosity down 
by a factor of 2, to log(L/L©) = —2.5. Both effects 
combined would be more than enough to solve the prob¬ 
lem. Other possibiliti es could include an unusually large 
radius, for example. Aller et al. (2013) also observed 
and reported the 41610-1913 B over-luminosity prob¬ 
lem, and concluded that J1610-1913 B does not look 
coeval with its host star. They also concluded the same 
for the five companions they observed in USco. In our 
study, however, we observed an over-luminosity only for 
J1610-1913 B, the other companions luminosities being 
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Figure 7. Spectrum of J1609-2105 b as observed with NIFS in the J band | |Lafreniere et al.1|2010[ | and NIRI in the H and K bands 
(|Lafreniere et aL||2008|l. The observed spectrum is compared with synthetic spectra trom tne tsI-bETTL models at logg of 4.0 (red) and 
o.U (^bluej and the JJrift-Phoenix models at logg of 4.0 (yellow) and 5.0 (green), for Teff=1600K (top row), Teff=1800K (middle row), 
Teff=2000K (bottom row). 
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Figure 8. Luminosity for different masses (labeled in units of 
Mini as a f unction of aee from the evolu tion models of [Burrows] 
jet al.|(|1997|l (red) andjAllard et al.|||2013|l (blue). The points with 
error bars mark the estimated luminosities of J1609-2105 b (pink), 
G06214 B (cyan), HIP 78530 B (green), J1610-1913 B (purple), 
and J1610-1913 Ab (grey). 


consistent with the 5-10 Myr isochrones in a luminosity- 
effective temperature diagram. If the luminosity prob¬ 
lem mentioned above for this companion were to be re¬ 
solved, the difference between its mass determined at 5 
and 10 Myr would be largely reduced. 

On the other hand, the masses can be estimated di¬ 
rectly from the evolution models by using the tempera¬ 
tures evaluated from the model atmosphere fits presented 
above. Specifically, we can find the masses for which the 
evolution models predict these temperatures. With this 
method, for an age of 5 Myr we obtain 0.007 ± 0.001 Mq 
for J1609-2105 b, 0.015 ± 0.003 Mq for G06214 B, 
0.020 ± 0.006 Mq for J1610-1913 B, 0.029 ± 0.012 M© 
for HIP 78530 B, and 0.15^!^;^! M© for J1610 Ab. For an 
age of 10 Myr, the corresponding masses are respectively, 


0.010 ± 0.001 M©, 0.016 ± 0.003 M©, 0.020 ± 0.005M©, 
0.030 ± 0.013M©, and 0.14lg;^jM©. The errors en¬ 
compass temperatures estimates from both sets of at¬ 
mosphere models considered here. The masses of J1609- 
2105 b and G06214 B are approximately the same when 
evaluated from luminosity or temper ature, and they 
agree well with the estimates m ade by jLafreniere et al. 


(2010) and Ireland et al. (2011) respectively. The mass 
of J1610-1913 B is significantly lower when estimated us¬ 
ing only its temperature, as was expected from the above 
comments. The difference between the two estimates ex¬ 
ceeds the quoted uncertainties, probably indicating that 
there is a real problem with the brightness of this object. 
For HIP 78530 B, this result is significantly higher than 
the previous value, but it comes with a large uncertainty 
and the two values can be reconciled. The mass of the 
close binary J16I0-1913 Ab seems to be in the stellar 
regime using both methods. 


5. DISGUSSION AND CONCLUDING REMARKS 

Our homogeneous comparison between the observed 
spectra of young substellar companions in USco and syn¬ 
thetic spectra from the BT-Settl and Drift-Phoenix 
models has revealed some interesting and systematic 
trends. First, the models do not succeed in reproduc¬ 


ing the details of the spectr a across the 1-2.4 ivm range 
simultaneously. As noted byjCushing et al.|(|2008|) in the 
case of early L-type dwarfs in the field, tne best fit in 
the individual bands typically occur for models of dif¬ 
ferent temperatures. At the temperatures providing the 
best fits in the J and K bands, the synthetic spectra 
have significantly steeper slopes in H than the observed 
spectra, both at the blue and red ends. An even more 
evident feature that is not reproduced by the models is 
the VO band in J. The VO absorption band at 1.06 /rm 
is only reproduced by the Drift-Phoenixs models, and 
generally only at a temperature slightly lower than the 
temperature leading to the best fit in other parts of the 
J band. The BT-Settl models simply fail to reproduce 
this VO feature at any reasonable range of T^s and log g. 
Similar conclusions about the VO feature and the fit of 
the H band for the BT-Settl models were reported by 


Allers & Liu (2013b). The alkali K I and Na I lines in 
the observed spectra are systematically stronger than in 
the models. Note however that we have not carefully in¬ 
vestigated the effect of metallicity on these features. The 
excellent agreement of the spectra of free floating BD in 
Upper Scorpius with our spectra (see Figure]^ provides 
yet another argument that the above trends are common 
features of young BDs and point to a real shortcoming 
of the models. 

Also, the best temperature estimates obtained by 
matching the broadband magnitudes and colors of the 
objects to synthetic magnitudes from the models are sys¬ 
tematically lower, by ~ 200 K on average, than the tem¬ 
peratures obtained from band-by-band comparisons of 
the spectra with models. This is also the case for the 
spectral fit applied globally (simultaneously across the 
near-infrared range), as it too is affected by the broad 
band colors. 

Beyond their use for testing atmosphere and evolution 
models, the wide low-mass substellar companions studied 
here are of high interest for the study of planet and star 
formation. In principle, such low mass companions could 
form like stars, through the collapse and fragmentation 
of a pre-stellar core, or as planets within a circumstellar 
disk, but their combination of low mass and wide sepa¬ 
ration poses a challenge to both processes. A formation 
in-situ within the circumstellar disk of the primary would 
require an unusually large disk, but a formation within 
a disk closer to the primary followed by outward mi¬ 
gration (from dynamical interactions) would be possible. 
The low-mass substellar companions studied here have a 
mass representing only 0.75-4% of the mass of their pri¬ 
maries. If these companions actually formed like stars, 
then it would imply that the fragmentation process can 
produce objects having only about 1% of the mass of the 
primary star. In any case, further observations of these 
systems using high-contrast imaging techniques and ra¬ 
dial velocity to search for additional companions would 
be useful to help understand their origin. For example, 
if they formed within a disk and were ejected outward, 
then a more massive object would likely reside in the 
system at a much smaller separation. 
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